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Abstract. The event-by-event analysis of multiparticle production in high energy hadron and nuclei colli- 
sions can be performed using the discrete wavelet transformation. The ring-like and jet-like structures in 
two-dimensional angular histograms are well extracted by wavelet analysis. For the first time the method 
is applied to the jet-like events with background simulated by event generators, which are developed to 
describe nucleus-nucleus collisions at LHC energies. The jet positions are located quite well by the discrete 
wavelet transformation of angular particle distribution even in presence of strong background. 



PACS. 21.65.Qr - 25.75.Gz - 02.30.-f 



Introduction 

Physics of high energy nucleus- nucleus collisions is much 
more complex than a mere independent superposition of 
nucleon-nucleon collisions. Already now the first results 
from RHIC on AuAu collisions at y/s = 200A GeV ev- 
idenced collective effects. The detailed analysis of RHIC 
results is presented in review papers of all four experi- 
ments BRAHMS [P, PHOBOS 0, STAR [3], PHENIX 
[4]. It is concluded that the matter, which is formed, is not 
described in terms of a color-neutral system of hadrons. 
This state of matter undergoes the most dense stage re- 
minding a liquid. It is called quark-gluon liquid or strong 
interacting QGP (sQGP). Manifestation of this new mat- 
ter can lead to such unexpected effects as multiple gluon 
minijet production, asymmetric and odd number of quark 
jets, collective effects. 

The most spectacular structure in the angular distri- 
bution of created particles are jets of hadrons well colli- 
mated along the jet axis. Comparison of jet production 
in AA and pp collisions or in the central and peripheral 
nucleus-nucleus collisions allows to estimate jet absorption 
in the partonic matter. This is one of the main signatures 
of QGP 0. 

Among suggested QGP signatures, the event-by-event 
topology plays a special role because fluctuations and cor- 
relations in angular particle distribution provide impor- 
tant information about the dynamics of a process. An ex- 
ample of such kind of structure is two particle angular 
correlations in the phase space at low transverse momenta 
Pt < 2 GeV/c which shows that jet-like structure exists 
even in soft hadron correlations and depends on collision 
centrality [3]. The collective azimuthal flow is seen in in- 



dividual events as their elliptic (" cucumber- like" ) shape. 
The ring- like structure of events has been observed [6,7, 
[5]. This phenomenon is interpreted as Cherenkov gluons 
PlOUllj or Mach waves [TW^ril] induced by a parton 
in a medium. 

The emission of mini-jets can be another origin of in- 
homogeneities in particle distribution. The clusters of par- 
ticles with relative small transverse energy are interpreted 
as mini-jets |15j . Mini-jet physics is the new trend that al- 
lows to investigate hadron fluctuations at low transverse 
momenta pt- At high px the jet production holds mem- 
ory about early hard partonic scatterings and jet mod- 
ifications due to the medium. Mini-jets at lower px are 
expected to have shorter mean free path in the medium 
and thus are more likely to dissipate, erasing initial cor- 
relations and changing substantionally their characteris- 
tics in a medium. Comparison of mini-jet event topology 
in central and peripheral nucleus-nucleus collisions allows 
to define more precisely the density and size of partonic 
medium and control a degree of its thermodynamic equi- 
librium [16j . 

The analysis of event structure by the wavelet trans- 
formation is very fruitful for such processes. The wavelet 
transformation was applied both to one- [TTUTB] and two- 
dimensional [19, 7, 20 particle distributions in nucleus-nuc- 
leus collisions. The discrete wavelet transformation (DWT) 
[7] and continuous wavelet transformation [19U20) were 
used. It was shown that the wavelet analysis allows to 
reveal fluctuation patterns. The "texture" of some events 
in AA collisions was investigated with the help of DWT 
[2lp2TfT6] . 

In the future experiments at LHC energies the jet topol- 
ogy of events can give the information about characteris- 
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tics of sQGP, produced in nucleus-nucleus collisions. It is 
necessary to reconstruct energy and position of jets with 
good resolution. Jet finding algorithms of cone type se- 
lect groups of particles within a cone of the certain radius 
Rj e t = \J (A77) 2 + (A(^>) 2 in pseudorapidity-azimuthal an- 
gle (77 x if) space (23j . This algorithm can't extract more 
complex structures in two-dimensional angular distribu- 
tions like rings. We show that this can be done by the 
wavelet analysis of a single event. 

Jet reconstruction by usual methods becomes difficult 
in central PbPb collisions at LHC energy [33] because 
of large background with (dN c h,/dy) y= o — 3000 ~ 5000 
charged particles. Therefore usual cone type algorithms 
are modified in different ways to take into account this 
background. In [25] the possibility of jet recognition in 
central PbPb collision in CMS detector with the help of 
modified UAl-algorithm was studied. It is shown that for 
jets with Ex > 50 GeV the recognition efficiency is close 
to 100 %. 

The two dimensional angular distribution in this work 
is simulated for central PbPb collisions at LHC energy by 
PYTHIA [26] and HYDJET event generators [24], which 
are aimed on the comparison with future LHC data on 
nucleus-nucleus collisions. 

With this goal in mind we try to find if the wavelet 
transformation helps to answer the following questions: 

1. Is it possible to distinguish different structures in 
the same event analysing the shape of their angular dis- 
tributions? 

2. How to find and separate the hadron jets with dif- 
ferent widths? 

3. What is a criterion of jet recognition over a huge 
background in heavy ion multiparticle event? 

We show that the method of the wavelet analysis al- 
lows to distinguish jets from rings and to find the jet po- 
sitions among a huge background even at rather low jet 
energies. 



The discrete wavelet transformation 

In many papers on wavelets 27 ; 28 , 29 it is shown that 
for expansion of a two-dimensional function f(x,y) in a 
wavelet basis one can use the one-dimensional expansion 
coefficients. Therefore here we explain briefly the basic 



g[n] = (— l) 1 n h[l — n] to the wavelet function ip(x). Then 
the orthonormal basis is constructed as a set of functions: 



{Vi,«(aO} = 2^ 2 ^x-n)} j 



(2) 



Wavelet function tpj^ n (x) expands 2 3 times with the 
scale changing and shifts along x axis by n. This allows to 
find features of f{x) not only on the whole x-axis (as in 
Fourier transformation), but also in narrow regions cor- 
responding to smaller scales. The restoration is complete 
due to completeness and orthonormality of the wavelet 
basis. Further we follow C. Mallat [27] . 

A great advantage of the discrete wavelet analysis is 
the opportunity to proceed with the so called fast wavelet 
transformation (FWT) using the coefficients h[n) and g[n]. 
Computing is done by the iterative procedure and is there- 
fore fast. 

Values of the function f(x) on the smallest scale are 
linked to the coefficients aj^ n = Oj[n], and further algo- 
rithm of transition to a larger scale is [27] : 



a 3+i[p\= E h[n - 2p\a j [n], 
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while the algorithm of restoration is: 
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(4) 

The sums in ([3]) and (J4| are finite so the finite number of 
coefficients h[n] and g[n] is used. 

In the two-dimensional case a separable form of wavelet 
functions is used |27j . Then orthonormal basis is con- 
structed as a set of functions {r£ nun3 , ^X ni> n 2 , r Em,n,} : 

,(xi,x 2 ) = 2 : >y(2?x 1 - m)^{2 : >x2 - n 2 ); 
,Jx u x i )=2?il){2?xx-nx)<p{2?x i -n 2 )\ (5) 
^Vn.naOi'^) = 2^(2^! - m)V(2^ 2 - n 2 ). 

Separable two-dimensional convolution can be defined as a 
product of one-dimensional convolutions on the rows and 
columns. Symbolically the decomposition and restoration 
of a two-dimensional function f(x 7 y) is shown in fig.[T](see 
27]). Here h[n] = h[— n]. Decomposition (fig. QJi) is done 
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concepts for expansion and restoration of the one-dimensiona|i rs t along strings of table a, (n , x 2 ), then along columns 



function f{x) 

There is the complete orthonormal basis of scaling 
{<Pj,n(%)} and wavelet {ipj, n {x)} functions with compact 
support in which one can expand any measurable function 
as a series: 
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One may choose wavelets which arc defined by finite num- 
ber of real coefficients h\n], called a filter. Coefficients h[n] 



Restoration (fig. [TJa) is done in the opposite direction. 



Wavelet analysis of complex structures in two- 
dimensional distributions 

The wavelet analysis is often called a mathematical mi- 
croscope since it allows to examine a signal with different 
resolution and to reveal structures on different scales. We 
demonstrate it here with some examples using Daubechies 
wavelet D 8 [2H]- DWT with Daubechies D 8 was tested by 



correspond to the scaling function ip(x), and coefficients the example of the one-dimensional function in [30] , 
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Fig. 1. Symbolic schemes of the two-dimensional DWT: a) decomposition and b) restoration 



In figure (figHK) the two-dimensional histogram of an- 
gular distribution in r\ x <p space is simulated from N = 
100000 entries, including the jet and ring-like shape as a 
toy example. The histogram binning is 128 x 128. The val- 
ues in each bin of the histogram determine the coefficients 
ctj=o[ni, n%\ for the smallest scale j = 0. The whole inter- 
val of the histogram along both axes 77 and ip corresponds 
to the scale j = 7. 

With the help of the wavelet analysis we try to re- 
veal different forms of irregularities in the histogram. For 
this purpose we use the wavelet decomposition of the his- 
togram and then restore its structure at different scales. If 
during restoration we set to zero the coefficients dj [71 1 , 712] = 
for large scales j = 3, 4, 5, 6, 7, then only a narrow peak 
is restored (fig. [2b)- If during restoration we set to zero the 
coefficients dj* [711,712] = for small scales j = 1,2,3 then 
the ring-like structure, which initially is present in the dis- 
tribution, is revealed and the narrow peak disappears (fig. 

lib). 

Now we consider another example of the function con- 
taining two peaks of different shapes: 



V(V - Vo) 2 + {<P - ^o) 2 
+ Aexp { -0.5 ( ^^ )2+ 2 ( ^^ )2 }. (6) 

The chosen parameters are a = 0.4 and A = 40 (the 
particular values for peak positions are irrelevant and cho- 
sen for peaks not to overlap). The peaks have different 
widths. Therefore they correspond to different scales in 
wavelet expansion. Just as in the previous example, we 
produce the histogram from this function (fig. [3^l). If we 
leave only coefficients dj [71 1 , 7x2] for j = 1 , 2 at small scales, 
we get the narrow peak (fig. [5b)- The wide peak is re- 
stored quite well if in the wavelet series we leave coeffi- 
cients dj [711,712] at large scales j = 4, 5, 6 (fig[3j;). 



Wavelet analysis of jet events 

Jets of hadrons are usually well collimated along the jet 

axis, so we search by wavelet method for narrow peaks in 

,2 

the distribution of the transverse momentum , P T . Jets 

depart 

with different transverse momenta px were simulated for 
pp collisions at LHC energy yfs = 5500 GeV by PYTHIA 
event generator |26j . Then the angular distribution of par- 
ticle transverse momentum was superimposed on the back- 
ground, simulated by HYDJET generator [53]. The back- 
ground is a sum of 100 events of PbPb collisions with 
multiplicity in the central rapidity region (dN/dy) y =o ~ 
4000. The jet production is switched off in the background 
events. 

Histograms are made by binning of the two-dimensional 
distribution within intervals \rj\ < 2.4 and \ip\ < it with 
the number of bins 128 for each variable. The sizes of bins 
are equal to Arj x A<p = 0.0375 x 0.049 , which are close to 
CMS detector granularities. The wavelet analysis is used 
to find the position of jets. The background is treated as 
the large-scale structure (in our case j = 2 — 7) . 

The following algorithm of jet recognition with the 
help of wavelet analysis is proposed: 

1. A given event distribution is decomposed by the 
wavelet discrete transformation. 

2. A smooth background at the large scale j = 7 is 
evaluated and then subtracted from the distribution. 

3. To remove the irregularities at the large scales 
the coefficients dj [711,712], j — 4 — 7 are set to zero. The 
range of small scales j < 4 corresponds approximately to 
the region yj (At;) 2 + (A^) 2 < Rj et — 0.3 from the center 
of the narrow peak. 

4. The coefficients at small scales (^[711,77.2], j = 
1, 2, 3, which are below a certain threshold, are set to zero 
in order to remove the sharp irregularities with small in- 
tensity. These cuts correspond to conditions dj [m , 712] < 
cutx [dj[ni, n2\)max- The parameter cut — 0.95 was found 
as optimum. 
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5. The event distribution is restored with the se- 
lected dj [711,713] coefficients. 

The event with three jets with total pr — 53 GeV, 
Pt = 45 GeV and pr = 22 GeV in their regions is shown 
in fig[4]before background was added. The sum with back- 
ground is shown in figO ■ The same event with jet recon- 
struction by the proposed method is shown in fig. [6] Here 
the smooth background has been removed by setting to 
zero dj \n\ , 712] at largest scale j = 7. The sharp peaks with 
small intensity in the background are removed by the cuts 
of small scale c£ 3 [ni, ^-coefficients below the threshold. 
Position of jets is well allocated on position of peaks in 
the spectrum of wavelet coefficients at different scales de- 
pending on the width of a jet. Fig. [6] shows that the peak 
intensities are decreased and their shapes are distorted. 

Only two jets are restored in this case. The third jet 
with px = 22 GeV/c is lost because it is compatible to 
background fluctuations. Another disadvantage of this se- 
lection is very strong distortion of the total transverse 
momentum values. The calculation of jet energy must 
be made by other methods (for example, by the cone 
method). But the jet positions are reconstructed quite well 
by DWT. 

The background under jet may be different and depend 
on jet position in 77 — ip plane. So it is not correct to sub- 
tract background, calculated as average background in r/ — 
ip space. Here we propose to estimate the background con- 
tribution under jet by calculating it inside a ring around 
the jet. Let we know a position of jet by wavelet analysis 
and select some region of jet as (Ai]) 2 + {Atp) 2 < Rj et 
. Then we can calcu late the total tr ansverse momentum 
in the ring Rj et < \J (A77) 2 + (A^) 2 < R r i ng - In our case 
we take Rring so that the number of bins under jet and in 
the ring are equal. 

Let's define the average transverse momentum < px >= 
IvKn J2k=i PT-k, where N bin is the number of bins in 
the region of jet or ring in angular distribution. The first 
(53 GeV/c) and second (45 GeV/c) jets have the aver- 
age < pr > equal to < pt >jet— 0.346 GeV/c and 
0.295 GeV/c. They are larger then those for background 
< Pt >back= 0.330 GeV/c and 0.246 GeV/c. For the third 
jet (22 GeV/c) in FigE] the average < p T >jet= 0.143 
GeV/c is less then for background (< pt >back= 0.209 
GeV/c). Thus the wavelet analysis can't find this jet. 

The best way to study jet events with large background 
by wavelet method is as follows. First, one finds positions 
of jets by suggested procedure. Then the background con- 
tribution is calculated in the ring region around each jet. 
The genuine energy of the jet is obtained by the subtrac- 
tion of the background contribution from the total energy 
in the jet region. 




Fig. 2. Ring-like and jet-like structures in angular particle dis- 
tribution: a) both structures, b) wavelet restoration without 
large scale coefficients, dj[ni, ni\ = at j = 3, 7, c) wavelet 
restoration without small scale coefficients, d 3 - [711, 712] = at 
j — 1, 3. View in the upper plane of figures is the distribu- 
tion projected on the plane. 



Conclusions 



Wavelet analysis allows to find jets and ring-like structures 
in individual high multiplicity events in nucleus-nucleus 
collisions. The two-dimensional discrete wavelet transfor- 
mation reveals well the irregularities with different shapes 
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c) 

Fig. 3. Narrow and wide jets in angular particle distribution 
a) both jets, b) wavelet restoration without large scale 
coefficients, dj[ni,ri2] = at j = 3, ...,7, c) wavelet restoration 
without small scale coefficients, dj[m, n2\ = at j = 1, 3. 



by corresponding choice of the wavelet scale j. Such struc- 
tures correspond to different effects in nucleus-nucleus col- 
lisions. The ring-like structure is impossible to reveal by 
well known jet algorithms, but this goal is achieved with 
the help of the wavelet analysis. The method is applied 
to the jet-like event with high multiplicity background for 
the first time. The DWT allows also to distinguish the 
bumps with different widths in two-dimensional angular 
distribution in the same event. It is very important for 
diagnostics of quark and qluon jets. 

Discrete wavelet transformation was tested on the par- 
ticle angular distributions, which are close to foreseen ex- 
perimental LHC data. Position of jets is well allocated on 
position of peaks in the spectrum of wavelet coefficients 




Fig. 4. The event with jet transverse momenta pr = 53 GeV, 
p T = 45 GeV and p T = 22 GeV, simulated by PYTHIA. 




Fig. 5. Sum of the event and the background, calculated by 
HYDJET. 




Fig. 6. The distribution of the event (see Fig(4} after restora- 
tion by suggested method of jet searching. 



at different scales depending on the width of a jet. Re- 
moving the wavelet coefficients at large scales allows to 
substract the smooth background. Removing the wavelet 
coefficients at small scales below certain threshold allows 
"to clean out" event from the sharp irregularities with 
small intensity. 

One must be careful with excluding the background 
by wavelet analysis. Its subtraction can decrease the jet 
intensity (energy, number of particles) and change the jet 
shape. It can be recommended for use only to find the 
jet position. The jet energy must be estimated by other 
methods. However, the jet position is defined quite well 
by the discrete wavelet transformation. 
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At the end we would like to stress that the wavelet 
analysis can be used more widely for studies of the gen- 
eral topology of individual high multiplicity events. Beside 
jets and rings, other structures of mini-bias events can be 
found as was noted already. Common features of corre- 
lations in particle locations within the three-dimensional 
phase space such as fractality and intermittency can be 
revealed 1511. 
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